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Abstract 

The angular part of the Schrodinger equation for a central potential is brought to the one - 
dimensional "Schrodinger form" where one has a kinetic energy plus potential energy terms. The 
resulting polar potential is seen to be a family of potentials characterized by the square of the 
magnetic quantum number m. It is demonstrated that this potential can be viewed as a confining 
potential that attempts to confine the particle to the xy-plane, with a strength that increases 
with increasing m. Linking the solutions of the equation to the conventional solutions of the 
angular equation, i.e. the associated Legendre functions , we show that the variation in the spatial 
distribution of the latter for different values of the orbital angular quantum number / can be viewed 
as being a result of " squeezing" with different strengths by the introduced " polar potential" . 

PACS numbers: 



* On sabbatical leave from the Department of Physics, University of Jordan 
^Electronic address: m ohammad@metu.edu. tr^ 
Electronic address: mustafa^purdue.edu 



1 



I. INTRODUCTION 



The Schrodinger equation for a central potential is a "classic " example in undergraduate 
and first year graduate quantum mechanics courses (see [IJ, for example). It is a good exercise 
in the separation of variables; a good example of the conservation of angular momentum, 
and a starting point for treating the Hydrogen atom problem. An important concept that 
is introduced along with this problem is that of the "effective radial potential", where the 
central potential that appears in the radial equation - after separation of variables-acquires 
and additional term that depends on the total angular momentum quantum number / (see 
below). One then has a one- dimensional Schrodinger equation with an effective potential 



The behavior of the radial wave function is interpreted in view of this effective potential. 
In particular, the new repulsive part, is interpreted as a centripetal barrier that attempts 
to throw the particle away from the origin, thus determining its shape for different values 
of /. The present paper analyzes the angular equation along similar lines (for the polar 
variable 6) and demonstrate that, viewed as a one-dimensional Schrodinger equation, we 
can identify an effective polar potential that shapes the angular part of the solution. The 
aim is to give the students of quantum mechanics a "feeling" of how the spherical harmonics, 
that constitute the angular part of the wave function for any radial potential, assume the 
distribution they have about the 2;-axis. We will identify a polar potential - more precisely, a 
family of polar potentials that depend on the magnetic quantum number m - that squeezes 
the particle towards the xy-p\ane with a strength that increases with m. The physical 
picture introduced, and the intuition employed, we believe, will also provide the students 
with a simple prototype of the intuitive thinking that is employed by practicing physicists 
in understanding more complex real-world systems. 

II. THE SOLUTIONS OF THE SCHRODINGER EQUATION FOR A CENTRAL 
POTENTIAL 

The Schrodinger equation 




(1) 



v^V' + vij = Eij 



2m 
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separates in spherical coordinates for a central potential by writing if) = R{r)Y {6 , (p) , and 
introducing the separation constant /(/ + 1) . The resulting radial and angular equations 
are, respectively: 



1 d f . ^dY\ 1 d'^Y ^^^^ . , 

sm 6 36 \ 86 J sm 6 ocj)^ 

The radial equation is brought to the Schrodinger form by substituting: 

u{r) = rR{r) 

so equation ^ becomes: 

+ { y{r) + \u = Eu (4) 



2m (ir^ |^ 2m 
with the term in the brackets being the effective potential we mentioned above, equation 

Now if we look at equation ([3]), the so called angular equation, we can do one more 
separation of variables by setting = G(^)$(0), then we will have the azimuthal 

angle equation. 



^1 + = (5) 



and the polar equation. 



1 d f , dQ\ rn? 



sin^-jTT = -/(/ + 1)6 (6) 



siYi6d6 \ d6 J sin^^ 
the solutions of ([s]) are 

$(0) = e'"^^ (7) 
and we require the solutions $(0) to be single valued, i.e 

<l>(0 + 27r) = $(0) 

which dictates that m = 0, ±1, ±2, .... It is worth to recall that the r.h.s of (|3]) is the total 
angular momentum hermitian operator represented in spherical coordinates. In general 
admits integer and half integer solutions for the eigenvalue /, which is the result of the 



angular momentum commutation relations. The restriction of / to integers is a consequence 
of restricting m to integers. 

Let's go back to equation (|6|, and write down its famous solutions which are the associated 
Legendre functions [Il[2](-/V/" is a normalization constant): 

0(9) =NrPr {cos 9) (8) 

where are given by 

I / \ 

prw = (i-:^')^(^^j Pii^) (9) 

where Pi{x) are the the Legendre polynomials which are given by the Rodrigues formula 

OH- 

The solutions \P[^{cos9)\'^ represent the probability of finding the particle at a certain 
angle 9, which means that they give the probability distribution of the particle about the 
z-axis. We will show below that we can view this distribution as being shaped by a polar 
potential. 



III. THE POLAR POTENTIAL 

Now we can bring the polar equation ^ to the Schrodinger form, just as we did with 
the radial equation. Defining: 

Q{9) = sin- ^9)y {9) (11) 
and multiplying by 1/2, equation ^ becomes 

where E = (l/2)(/(/ + 1) + 1/4). 

This equation can be thought of as a one dimensional Shrodinger equation (with 
h = m = 1 ) for a particle confined between and vr, and satisfying the boundary conditions 
y{0) = y{7T) = 0. The term in brackets is a one dimensional potential, more precisely 
it is a family of potentials depending on the choice of the value of m. These potentials 
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FIG. 1: Graphs of the potential for several values of integer m 

are special cases of the so called Poschl- Teller potential [31 H] . To gain a closer look at 
the potential, we plot it for various m-values in Figjl} As the figure intuitively suggests, 
the potential, which is infinite at the boundaries - thus banning the particle from going 
beyond the boundary- attempts to confine the particle about the 2;-axis, the higher m is 
the stronger is the confining strength. For the lower m values the potentials are very close 
to an infinite well with width tt, as we increase m the potential squeezes the states resulting 
in more confinement. 



Let's take a look at the solutions now. The solutions for this family of potentials can 
be worked out in several ways, a nice algebraic way is to use Infeld and Hull factorization 
method [3]. This family of potentials enjoys shape invariance symmetry, which can be 
exploited to workout the eigenvalues algebraically. Supersymmetric quantum mechanics 
techniques can also be used to work out the eigenstates [6| . We will not use these techniques 
here, rather we will work out the solutions from the solutions of the original angular equation 
(|6|. If we restrict the values of m to nonnegative integers, the solutions, upon employing 



the transformation (11) can be written directly from the solutions of the original equation 
and they read 

CW = iV-sini^P„7^(cos^) (13) 

where -P^^ ^^^5 again, the associate Legendre functions. These solutions are for non- 
negative m, those for negative m are - as is well-known [U [2]- directly proportional to these 
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FIG. 2: I and m lattice 

solutions. So, without loss of generality, we will restrict ourselves to nonnegative m values. 
In the original equation ([6]), the solutions are conventionally labeled by a pair of numbers 
(m,/) and subject to the condition m < I. The pairs {m,l) constitute a lattice as shown in 
Fig|2] In the conventional interpretation of the solutions of ^ the number I is associated 
with the angular momentum, so the usual way to look at the grid of solutions is to fix the 
angular momentum by fixing / and then changing m, so all the solutions with different m's 
and same / are degenerate and correspond to the same angular momentum. This amounts 
to moving right (or left) on the grid. But what we are doing here is different, we are fixing 
m and looking at the solutions of a one- dimensional potential at a fixed m value. When 
we fix the value of m to a specific integer we get all the solutions for that one dimensional 
potential. As one can see, the first eigenstate for m = 1 is at / = 1, but it still corresponds 
to the ground state of m = 1 potential, which will be labeled as n = 0, where n is the 
quantum number which labels the states of a fixed m potential. The same thing happens 
with all other m-values, thus we are moving up the grid shown in Fig|2} It is always the 
case that the ground state of a fixed m potential coincides with a / = m, therefore, the 
relation between / and n is / = m + n. Consequently, the eigenergies can now be written 



down directly from those of equation (12); they read: 

K = lin + m + ^f (14) 

Here, higher values of n, i.e. excited states correspond to higher values of the kinetic 
energy of the confined particle. That this corresponds to higher /-values -in the conventional 
view- is natural, as higher kinetic energies mean higher angular momenta. 
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Now, we show how to relate the behavior of the potential to that of the eigenenergies, 
eigenstates and the original associate Legendre functions. Let's look again at the graphs 
of the potential for several m in Fig.Q. Here we note that as m increases the bottom of 
the potential is shifted up, and this is responsible for increasing the ground state energy 
with increasing m. The second observation here is that as one goes for higher values of m, 
the potential tends to confine the particle closer to the mid point of the potential = f), 
which means confining the particle closer to the xy-plane. This is confirmed by looking at 
the behavior of the ground states for several m. In Fig. ([s]) we plot the ground states for 
various values of m just below the corresponding potential. It is clear that as m increases the 
ground state gets squeezed closer to middle of the box, giving rise to more confinement of 
the particle. The last line in Fig. ([s]) shows the corresponding associated Legendre functions 



which are related to the eigenfunctions through the transformation (11). It is quite obvious 
that the particle is more likely to be found near the xy— plane as m increases. Since higher 
m means higher angular momentum, this squeezing of the states towards the xy— plane can 
be explained semi-classically in terms of the increase of the moment of inertia of the particle 
about the 2;— axis. Finally, while the examples we have considered above are for the ground 
states only, the same arguments apply equally well to the excited states. Fig. Q shows 
the first excited states; n = 1, for three m-values. It is clear that as one goes up in m, the 
particle gets more squeezed towards the xy-plane. 

In conclusion, we see that we can view the angular distribution in the wave functions 
of a particle in a central potential as being a result of confinement by an m— dependent 
"effective" polar potential. 
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FIG. 3: The first row shows the potential for three values of m, the second row shows the ground 
state solutions for those potentials, the last line shows the polar plots for the Legendre functions 
with the argument cos 6 for the corresponding values of m and / 



m=0 ,n = l 



m = 10 ,n = l 
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FIG. 4: The first excited states for m = 0, 10, 30. The excited states get more confined about 
6 = it/2 as m increases 
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